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ABSTRACT 

A detailed examination ox the literature pertaining to subccoled 
critical heat flux x^as conducted. Subcooling, system pressure, mass 
flow rate, system geometry and flow stability were determined to be 
the important parameters in determining critical heat flux. The effect 
of each of these is discussed in detail. 

Various schemes for correlating experimental critical heat flux 
were subjected to critical examination. Form, content and comparison 
with experimental data were the bases for evaluation. The correlation 
of Griffith was found to compare most closely with the experimental 
data over a broad range of variables. A linearized correlating equation 
was developed for the MIT low pressure (30 to 90 psia) data, including 
that determined as part of this study. 

Critical heat flux was determined experimentally for pressures 
from 30 to 90 psia, mass Hew retes of 3.69 x 10*° to llj.76 x 10° 
lbra/ft^-hr and diameters of 0.0lt7 to 0. 2l*2 inches ever the subcooled 
region. The experimental data obtained is in general agreement with 
that of other investigations. This low pressure data, when combined 
with the results of previous investigations, provides a fairly complete 
picture of subcooled critical heat flux for tubes up to 0.25 inch 
diameter. A graphical interpolation scheme is presented by which it 
is possible tc obtain predictions of critical heat flux which are 
generally superior to that of the mathematical correlations. 
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I. IWTRC XTTIO’-t 



In the continuing search for devices of hi; her energy flux _<3nsity, 
mere and more emphasis is being placed upon boiling as a means of 
increasing heat transfer rates with ’-•’oderate temperature iil feronres. 
Current developments in high energy electromagnets, high performance 
electron tubes and rocket motors require very high heat fluxes in 
cooling applications. Although heat fluxes well in excess of One 
Million Btu per square foot per hour can be obtained, heat fluxes 
cannot be increased without limit. For a given set of flew conditions, 
there is a limiting Critical Heat Flux (CHF), at which the transition 
from nucleate to film boiling occurs. At this point, the wall temperature 
increases dramatically, usually resulting in catastrophic thermo-mechanical 
failure of the heat transfer surface. One also finds references in the 
literature to this as "burnout”, and "the boiling crisis". 

Often, the most desirable mode of boiling for heat transfer 
applications is that of surface, or subcooled boiling, whe~e the average 
or bulk fluid temperature remains below that of saturation. In this 
case, vapor formed at the surface condenses when it comes in contact 
with the cooler main stream. The Pressurized Water Reactor uses water 
at high pressure (2000 psi) as coolant end moderator . The high pressure 
allows sufficient coolant temperatures for further heat transfer and 
power generation, while maintaining bulk temperatures below saturation. 

The maintenance of subcooling is important in reducing void fraction 
and retaining continuity of the moderator. The electromagnets and 
electron tubes are generally coded by low pressure systems where the 
expense and complications of high svstem pressures can be avoided. 
Utilization of subcooled boiling eliminates the complexities of two- 
phase flow, as well as the mechanical complication cf condensing 
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equipment ; however, the* CIS- problem is rti.ll p -sent in r decoded 
boiling. 

Cvor the poet several years, numerous studies have been concerred 
with this CHF problem. A great deal of experimental data is available, 
under various conditions of flow. Invest!, voters have usually selected 
certain parameters and investigated the effect of varvinp these within 
certain limits. The majority o< the previous experimentation has 
been conducted with uniform heat flux and round tubes. This investigation 
has been limited to this region. 

The first phase of this investigate on consisted cf a search of the 
literature to determine the flow conditions which had been investigated. 
Techniques, generally in the form of correlating equations, have been 
proposed for predicting CHF, by almost evervene who has conducted 
investigation in this field. One phase of this study has been concerned 
with a critical examination cf the techniques and their form and 
agreement with available experimental data. 

Another phase of this investigation consisted of completion of 
some previous experimental work. Experiments were conducted at 
selected tube diameters, at low pressures and various flew rates. 

Data thus obtained, when combined with that in the literature, provides 
a picture cf CHF, for tubes of small diameter and flow rates of interest, 
over the complete range of pres-ures, from 30 osia tc just below critical 
pressure. 
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This iis a two part stud"; c su.-'.ey c' t. 5 (i -> 1 ,r ~e 

and correlating equations ar.c experiments < ict. r»an>- -ic c. "H. •!c"_ 

Krat Flux ( CHF ) under certain c nd: t" cns. 

First, citations oi all oi the current literature cmre^ne.' V:‘th 
Subccoled Critiftal Heat Flux were obtained. From these, the majority 
of the literature was reviewed, and experimental data err ore csed 
correlating equations were reccriei for further investigation. i 
largo part of the work in this field has been accomplished in Russia, 
consequently certain information is available only in the Russian 
Journals, untranslated. It was usually possible to extract data and 
proposed correlating equations without complete translation of the 
particular paper. 

A listing of proposed correlating equations was prepared. These 
were compared with actual experimental data from various sources, 
including the experimental phase of this study, in an attest to 
determine the validity of the correlation schemes. 

Experiments were conducted to determine CEr for predetermined 
values of the important variables. The dependence of GET uren 
subcooling was determined for three fixed values of mass flew rate, 
with c cns cant system pressure, tube length and diameter. This data 
completes the low pressure end of the study of Ornatski and Kichigin 
(b, f> and 6). Then, using one oi the same flow rates and pressures, 
and holding exit subcooling constant, dependence of CHF upon tube diameter, 
and upon tube length-diameter ratio was determined, Next, the variation 
of CKF with changing entrance forms, with the ether variables held 
constant, was investigated. Finally, tube diameter was changed to obtain 
comparison data, and to determine the maximum CHF available at the 
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larger diameter with the experimental Apparatus. 

All of the experimental data relating tc Subcooled Critical • -at 
Flux that has been obtained within recent years at M IT was collected, 
converted to similar format and dared on machine data cards in 
anticipation of future studies. 

A detailed ^-description of the experimental apparatus and procedures 
employed is contained in Appendix: C and E. 
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'ArTCPP ~ ,r rT.vP fr T''C CRITIC A.T F frLUX 

Although there has been a rest deal cl research and i vestige tion 
in this field, particular? in the Ir-s- decade., there exists tods r no 
clearly dellnci, universal, mathematical dencripti.cn of the process 
which results irV a limiting CKF. Pool Belling is itself a complex 
problem, the additional conditions of ^ubccoled, forred flew tending 
to cloud the picture even rore. Attempts at dimensional snail’s is, 
and treatment of the problem from basic equations have, in the rain, 
produced little information, usually yielding long, complicated 
functional relationships and involving simplifying assumptions net 
necessarily borne out by experiment. 

It is generally accepted, however, that subcooling, system 
pressure, mss flew rate, flow geometry and stability arc the major 
factors in determining CHF. These variables evidence various decrees 
of interdependence, and in discussing the effects of each, this 
interdependence must be berns in rind. 

Effort of Subccolin- 

Subcooling is essentially a measure of the amount of pre-heating 
required to obtain boiling. This is generally measured as the 
difference between outlet bulk temperature, or mere precisely enthalpy, 
and that at saturated renditions. This is reasonable for uniform 
heating, in that failure normally occurs quite near the downstream 
end cf the heated section. From another viewpoint, subcccling is a 
measure of the cooling effect cf the main stream. Thus CRT increases 
■with subcccling. There is a tendency in the literature to treat CHF 
as linear v/ith subcooling. This is not the cose, but it is in general, 
not a grossly inaccurate approximation. In the region of low values of 
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cubcooling the data air a ' st ‘ -r- < i- ^**rr „v, t a nini* ral 

CIT is -r c-'l.e. . 

f J fe r t r f r- 'C tr r , ■ lS. uro 

System ircsrurc pla^s at it “ ortan« role . i .e^cr - .i. , CT1 , altfecjp/- 
tic effect is les: era ratio than in Feel hcilinr . In a reccrt .‘.J c-nt 
study, it is stated that a variation of ays ,e: treats c iron iO psia 
tc 10^0 psia has little significant ef: ect on CM? this in tubes ranging 
from. 0.3 inch to 0,8? inch in diameter. However, this is net borne 
out by the other expo imental data ( 9 ). A pressure increase iron 
600 psia to 1000 psia results in a reel action o i CM ci as much as 2%. 
Further increases in pressure above 1000 psia result in further 
reductions in CHT. At lower pressures there is disagreement ever the 
effect of pres cure, McAdams (32) and Gunther (33) report no effect up 
to 163 psia, while Marshak, et al. (31) shew a pressure dependence over 
the entire pressure range considered. Eerg’ es (ll), reports a 15£ 
increase in CHF , for an increase in pressure frer 20 to 85 psia. 

Chirkin and Iukin (21?) demonstrate that increasing pressure, in the 
lew pressure region, increases CHF. Their pressure effort is nuch 
stronger than that cl other investigators, however, it is non-: felt that 
most oi their data suffered fror instability problem which could have 
clouded the real picture. 

The experimental data decs suggest, however, that ^ or fixed geonetrv, 
flew rate and subcocling under conditions of stable flow, there is sene 
pressure for which CHF attains a maximum value. The pressure at uhich 
this maximum value cf CHF occurs is net clearly defined, being dependent 
upon mass flew rate, and possibly the other variables as well. It is safe 
to say theu h, thao under conditions of reasonable lieu rate, this maximum 
occurs somewhere between IiOO and 600 psia. Unfortunately data in this 
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high values of sxibcocling (4h sa ^. greater th»n 200 TJ/lb), CIT becomes 
somewhat independent cf subccoUng. ' thor high pressure lata shows 
this sane general trend (l, $ ard (.), This w<uld, of crura*, reouira 
further investigation to be throughly substantiate ’. 

Effect cf M ass I lew Tate 

Generally, in the subcode! region, increasing mas- flow rate 
results in an 5 a* reused Oil. This is easily visualized since the 
increased flow sweeps the bubble away frer tie boat transfer surface 
ncro rap*' d'y, replacing then wiuh relatively cold waver. however, the 
rate of increase in CIT due tc increased flow rate is pressure dependant. 

At high pressures, the effect cf increased mass flc. rate becomes less 
pronounced. 

Effort c f Geometry 

Tho majority of the research end study i, this field has been 
1' rated to round tubes and annular fluid pas r. a gee. The gross erects 
of the other variables however, c re much the ear' i o~ a ay cross section (8). 
In general, for iven flow rate, pressure and subcceling, CVi is greatest 
for c circular tube ard successively lower as ohe- geometry becomes core 
complex. "Swirl flow" is one exception to this. The ratio ci heated 
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perimeter to otal vetted r*5-’inote~ h-'f bton rur/ertr 
the 03se of "swirl flew" (1?). 

Investigators disagree whet r or there is a ran where meter h's 
no effect, Pcroshchuck and Lar-tsman in a recent study (9), conducted 
tests for fccr tube diameters over rc.s cf the high pres'-ur" ran; e , 

(735 psia to 2500 psia), with exit conditions from high cubrccling 
(^^sat^^OO BTU/lbjjj) to bulk boiling (x exi ^-0.2). This study showed 
CHF decreasing consistently with increasing channel diameter. Earlier 
investigations show decreasing CIIF with decreasing diameter. Beth 
Chirloin sc lukin (2k) and Kafenhaus and Bocharov (26) report this. 
Again, this is now felt to have been in a large part due to flow 
instabilities in their systems. ^uFont (8), Pirshak (3) and 
Zenkevich (11* ) note no diameter effect. Bergles (Ll) notes a strong 
inverse dependence cf CKF on diameter for small diameter tubes, but 
he shews the Zenkevich Prediction (lit) as a possible asymptote for 
diameters greater than 0.3 inches. 

With flow rates high enough to eliminate the effects of natural 
convection within the passage, channel orientation has little effect 
on CKF. 

The effect of heated length has not been completely determined. 

In general, the length, or length-diameter ratio, is not important if 

the length is long, or L/D is large. The popular conception is that 

CHF is determined by purely local phenomena, thus as long as the flow 

is fully developed and entrance effects play no part, the length is 

not significant. But the question of how long is long enough has not 

been satisfactorily resolved, partly because of the interact icn of the 

variables. In general, the Russian studies consider an L/D ratio of 

10 to be large enough tc eliminate the length effect. This conclusion 
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has not been ccrp? etly verified. Perglea (11 V, ■' n » Ip*-? Tt»esu»r study, 
shewed an l/ r dopG n dr-nre existing -Per T./ r :n th* order c" I''', Pcrcshr’r ucl. 
and Fried (l8), r resent data for tubes of 0.1? inch di* meter at IvOO osie. 
that she;,"- nc dif r ercrce in CHf for I./ r of 1L , f 0 and 100, over a fairly 
braod rar.-e cf exit renditions. "MPont (8) suggests consider?- on of 
the to-n^oraf ure 'srro^ient ' cv,n ' he horted tube. The tempo " urn ! ij .ore nos 
is of o curse dictated by the First Law of Thermodvnanics , ■5 /s h»Iv.u this 
gradient really reflects is a length-diameter consideration. 

In general, CHI decreases for increasing I/P, but for high pressures, 
the dependence of C!T upon L/D becomes less significant for values of 
L/ n greater than 30. 

Effect of new Stability 

This refers to pulsations in flew and Pressure, resulting from 
compressibility effects in th-» heated sections. If the’-e is an expander 
or accumulator ups team, cf the heated section, there is a distinct 
possibility of producing a flow oscillation which will nericdically 
reduce the flow ne»r the tube cutlet and p r •*-bops cause a premature 
failure of the heat transfer surface. Bermles (11) points cut that there 
are two tvpes of instability* flow oscillation due to r empress ibility 
and flow excursion due to system characteristics. Flew conditions under 
which these pulsation* do effect CM are considered unstable. The result 
of unstable flow is a drastic reduction in *~HF . vith unstable flew 
resulting i* heat fluxes as low ns ?.Q* of the stable value. Otrble flow 
usually can be insured for subcooled boiling fcr installing some sort of 
throttling device (orifice, valve, etc.'* in the cistern just up-strean of 
the heated section. Doroshchuck and Lents'-*" used a pressure drop of 
lt70 psi in their high pressure studies (9, 10). Perries (11 ) reported a 
pressure drop of approximately 18? psi in his lew pressure work. 
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Further d^ soussiors o: the eta .r>. t' r pre * •■ *’ rr r i. i i 
Kef, 7, 11 >.r’ ?8, I'he stability Tf'lr i« rartlculcr-'’ ar ^ i n 

system roc o-*ini»'g ho it axel angers »nth parallel jbes from '’c^cq 
inlet and exit headers, where s subs on; ti *’1 amount c? €■ r< .rr; »1 -le 
velar© ray exist. In any e/eet, toe inr all- 1 - * on of • n- assure urin> 
levies, and Dhe^remiltant less of down-stream svst.' , n prs na; , j n «*o re 
than iustl^ied by \,i ■■ increase in CKi , 

Effect of Secondary factors 

Several other 1 actors have been investigated tc deter ine how they 
effect CH: , 

Surface rcaghnes- , for reasonably smooth (corarcrcial quality) 
channels has little or no eiiect (7, 13). 

Dissolved pas in the flow stream, in sr-ell amounts , has a negligibly 
6 -"11 eft ere on hr' (17). 

For heat flux unevenly distributed around the heated -oorimeter, 
peak heat flux was found tc increase. For heat flux unevenly distributed 
aloof the length of the heated channel, qualitatively similar results 
have been obtained (7). Tests have been conducted under conditions whore 
the boat flux was forced into a cos'ne distribution along the tube and 
CHF was shown to be reduced by seme 30% (11), 
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General ^iscus -ion 

t'omercus schemes have oeen cro^osed as ’orrc • a , ~ r . 

bv which rT v can supposedly be predicted. The earlies cf the •’"© 
of a never function fern: 

* (q/A) cr - A* &c ^ *£ 

whore 3, n and I are the important variables, cr combinations of them, 
and c, e and g are expo rireo nt a 11 v . eter~i ned exponents • usually 
eonreloti cne cf this + eT’ tc neglet l.h in '■’’artier bet :aen 

variables; in cth^r c-'s.-^f the experent? t barrel von *> ” o unctions of 
the controlling pc '■•'■meters . Another bread caterer/ Lac lueses the 
"superp sit 'erf 1 correlations, thos'* r- reining pool boiling data with 
forced cc nvcction heat trailer information. There correlations are 
often Iona one difficult tc use. 

Several attempts have boon made at analytical treat rent , Zenkevich 
(lb) and Fyabov and 1 eraine ( 36 ) with dirersioral analysis and Griffith 
(13, 1?) with a theoretical rral^ria. Other eqiictic v r were developed 
by "brute frree" data analysis, '$ A r>~ c nr put i n, ycMnes. ^aPcnt 
f avers a ferr: 

(q/A) cr « c (1 ♦ r x .A) (1 * r 2 .s) (1 = r 3 ,^) 

where A, B and I) are the influencing variables, (pressure, subccclinp, 
diameter, velocity,) and C, C-^, Cg, etc. a^e constants. 

T n general, the ccrrelnticns are limited tc: specific ranges of the 
variables, departing radically from the experiment?! data in the excluded 
areas. Tn sore o-ses, c ,r, e f i nds <'crr^laticns developed fro* a p*e*t of 
an exmrireentri study, using cert-i 1 ' d~ip on-: excluding ether. This 
technique is usually accompanied bv a statement to the of tor t th?t. so 
rarry percent cf the experimso ;al results agrees withi ■ plus cr ^-inus so reach. 
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In cases such as this, the potential v *r is obi' ed to examine the 
rj." r>rore r ’ure as closely as he does the correlatn.cn itselx. 

As a re sult ci this - tate of affairs, one phase of this study 
was concerned with examination of correlating equations. The listing 
of correlations is not all inclusive. Several correlations were examined 
and discarded, generally due to their talked similarity tc others. 

Of the correlation schemes examined, the relieving are those 
which agreed reasonably well with the experimental data, or illustrate 
some point. 

Details of Particular Equations 
Bell; (ref: 21) 

(q/A) cr » 0.1*27 x 10^ x 

Range of values; (q/A) cr * Btu/hr-ft^j p » 2000 psiaj G ® .2 x 10^ - 
5.0 x 10^ lbm/hr-f ; ^hg a ^ = 0 - 130 Btu/lbm 
This correlation is compared with experimental data (figure 9) 
for conditions very nearly identical, from three different sources. 

In obtaining this data from the various references, a certain variation 
was unavoidable j two people do not conduct tests at precisely the same 
conditions. However, the valuation was kept to less than the expected 
experimental error. The data from reference 15 was obtained in 195u, 
prior to the general acceptance of the need for stability, and may be 
of questionable value. 

It can be seen that the correlation follows the general trend of the 
data, with variations of 15? to 30%. 

Jen3 & Lottes t (ref: 17) 

This correlation was based on data taken at UCLA and Purdue over a 




range of pressures from 500 - 3000 psia. The values of K and n were 
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determined on the basis of the data such to make all deviations positive. 
(q/A) cr - K (C/I0 6 ) n *(^T sa ^ * 22 

UCLA Data correlated tb ~0$, *23%» D ” .26 in. 



_2_ 


K 


n 


5oo 


.617 


.16 


1000 


.626 


.275 


2000 


.Ui$ 


.50 



Purdue Data correlated to ~0% 9 +6C$, P ° ,H3 in. 



_L_ 


K 


n 


1.000 


.915 


.275 


2000 


•5U5 


.50 


3000 


.30 


.725 



The correlation was plotted for data (figure 9) of approximately 
D ® .ICO in., G = 2.1 x 1(/ lbm/hr-ft c , at p • 2000 psia using the 
constants and exponents of the Purdue data. 

Phrabov and Berzina: (ref: 36) 

1) (q/A) CI . 58 (q/A) pb + 1*03 "sat* ^sat*^ (/ 1 ~/° 1*^ 

(j 1 £gpv)' J ^ <r • lj °5 g *3 

Range of values: D * 2 tm a .079 in. , £ ® 18 - 590 deg F, 

I/O - 20, p - 10 - 200 Atmospheres - Uj7 - 29i0 psla, G - 3.69 x ID 6 - 
22.11i x JO 6 lbn/hr-ft 2 

2) (q/A) cr - (q/A) ph * 257 V-^W 67 kg Vt‘ 33 (V 1 ' Pv) l' 3 



Range of values: D ^ h mm (.158 In.), T sab * 16 - 1*30 deg F, 

I/D = 10 - 220, V « 3 - U8 ft/sec, p - 70 - 210 atm «* 1030 - 3090 psia 
Both of the above correlating equations were developed from 
dimensicnles3 equations. In each case (q/A) ^ is the same and is obtained 



1 ? 



from a plot of pool boiling OFF derived as a function of pressure, 
correlating equations were derived from a dimensionless equation* 
1) K aq = 3710 K w * 6 K^x* 67 Pr“*^ (^ - 2rm = .079 in.) 
or 2) K Aq «• 2370 %* 6 K aT * 67 Pr“ #il (D * lias * .158 in.'s 
where K^ q • (4q/A) cr h f ^ J go 

1 3at ^°1 “(°v)s 



*W - V 




I* 



Kj* A sat h f g p v J gp 
g o" T sa t -p v 

It is assumed that (q/A) cr » (q/A)^ + ( Aq/A) er 

The 2 mm (.079 in.) correlation was derived from the data of Ornatski 
and Kichigin (b> 5» & 6) by the usual technique of letting the data 
determine the coefficients and exponents. Generally the data falls within 
- 30$ of the correlating line. The correlation was also computed for the 
2 mm data taken during the experimental part of this study. In this case 
the correlation did not follow the basic trends of the data and was in 
error by as much as 100$ at 3 one subcoolings. 

The correlation for D ^ b mm (.158 in.) was computed and plotted for 
G ■ 2.1 x 10^ lbm/hr-f , p « 2000 p3ia and compared with the data far 
Poroshchuk and Lantsran 1 s (9) b ram and 8 mm tubes. The correlation fit 
the 8 rro (.316 in.) data very closely , but was 15$ - 30$ low on the 
b mm data. It was not plotted on figure 9 to avoid confusion as it 
fell generally cn top of the correlation of Zenkevich and Subbotin. 
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Zenkevic h : ( r ef : lii ) 



(q/A)ct * kfg <rs 80 ^ 3 x tf* + X ^ h sat] x 3.0 ^ 

Range of values: p » ll* - 165 psia, Ahg a t * 20 - 282 Btu/lbra, 

G = .22 x 10 6 - 9.96 x 10 6 lbiq/hr-ft 2 , D * .16 - .50 in. 

% 

This correlation is linear with subcooling as co’nnonly noted for 
Russian analyses. Specific instructions as to the temperature to 
evaluate viscosity at are lacking. Using the fluid saturation temp- 
erature yields an equation linear in subcooling, while using fluid bulk 
temperature produces a slight curvature in the resulting plot. This 

•y 

correlation is plotted in figure 8 for both conditions and it can be seen 
that the equation does not really reflect the experimental results. 

Zenkevich and Subbotin: (ref: Ik) 

(qA)cr “ 397 G^ AT sat «33 (v g /vf g )” ‘ Btu/hr-ft 2 

Range of values: AT sa t = 18 - 180 deg F, D*l*-12mm sa .158 - .1*72 in. 

p » 11*0 - 210 atmospheres * 2060 - 3090 psia. 

This correlation was plotted for G * 2.1 x 10^ lbra/hr-ft2 for 2000 
psia and 2500 psia and compared with the data of Doroshchuck and Lantsman 
(9) for tubes of 1* ram (.158 in.) and 8 ram (.316 in.)® At 2000 psia the 
correlation does an almost remarkable job of predicting CKF for D * 8nm® 

At 2500 psia the correlation is very close to the experimental 8 ram data 
although approximately 2055 low at low subcoolinga. The correlation 
does not fit the 1* ram (.158 in.) data as well, with differences of 15$ - 
30$ at both pressures. Figure 9 shows the correlation at 2000 psia for 
the 1* ram diameter. 
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Gambill: (ref; 30) 



1 + /pA ~ c p A -sat | 

in 7 n?nrf£" 

+ K» (k/D) Bc n Pr” (t a - t b ) B0 

where K * .12 - .17, K’, n, h n depend on the non-bcd2±ng turbulent 
flosi correlation used. 

Range of values: V = .05 - 17h ft/sec, p * h.2 - 3000 psia, 

AT sat “ 0 - 506 deg F, a = 1 - 57,000 g's 

This correlation is of the superposition form where a subcooled 
pool bdling burnout term is added to a non-boiling forced convection 
term, typically the well known McAdams correlation for turbulent flow. 

Gambill leaves some latitude for the evaluation of the constant 
coeffiecients and exponents. Evaluation of the wall temperature (ty) 
in the forced convection term is done by an auxiliary plot of AT sa +, 
vs. T sab /T cr it based on Be math's correlation of wall temperature at 
burnout (16). This is explained in detail in the section on Bemath's 
correlation. The correlation was plotted under different conditions 
with tire following results: 

1) D * .2h2 in., G * 3.69 x 10^ lbm/ft 2 hr, p * 90 psia (figure 8) 
The correlation fell about 35$ lower than the experimental data. 

2) D = .079 in., G - 3.69 x 10 6 lbm/br-ft2, 7.38 x 10 6 lbm/hr-ft 2 , 
& lli .76 x 10^ lt*a/hr-ft2, p « 90 psia (figure h) ; The correlation 
fell a3 much as 100$ lower than the experimental data. 

3) D a .180 in., p * 2000 psia, G 3 2.1 x I06 lbm/hr-ft (figure 9) 
The correlation is conservative (lew), but does not show the trends 
of the data completely, and is in error by as much as 255. 



(q/A ) cr * K hfg ^r gp a ( f 1 v) 
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Be math: (ref; i6) 



(q/A) C r " ^bo “ ^b'bo 

where; h^ c * 5 710 + U8 V D”*° Btu/hr-ft^-deg i 
ate: t*)^ - 32 ♦ (?/5)^5? ln(p) - $U(^ - vj J*t S 



Range of values ; Round tubes, p * $00 - 3000 psia, 

51* ft/sec, where D is measured in feet, and p in psia. 

This correlation is based on a Uniting teat transfer coefficient 
and a correlation of wall temperatures at burnout. The correlation 
was plotted (figure 9) for G ■ 2.1 x 10 C lbn/hr-ft^ at 2000 psia, and 
D * .180 in. for comparison with experimental data. The correlation is 
slightly conservative predicting this data within approximately 5$» 
Bettis : (ref: 17) 



(q/ A ) cr B 10 6 x 0.28( hh/10^ )“ 2 *^ ( 1 + G/lO 7 ) 2 e'*' 00121 '/^ Btu/hr-ft 2 
Range of values: G * 0.2 x 10^ - 8.0 x 10^ lbm/hr-ft' , p = 1850 - 

21^0 psia, Ah^ - 0 - 180 Btu/lbn, l/D - 21 - 365. 

This correlation is compared with experimental data on figure 9. 

It can bo seen that it follows the general trend of experimental data, 
with variations of 15$ to 1*0$. 

Griffith: (ref: 13, 17) 



(q/A) cr 



f( P/p 0 ) (F) 



g < h g * h b ) 




.33 



(F) ■ 1 * (ft) b x XO" 6 ) - O.Xit Ki + 0.5 x 10" J (Re,, K, )i 

** ’ V-A * % - A e , x (fct - *b > . 

yUl /*g h fg 

and, f( p/p c ) is presented graphically in the reference. 
Range of values: p * Hu 7 - 3000 psia, V * 0 - 110 ft/sec, 

Z>T sat - 0 - 280 deg F, 



21 















■ . , 

i x J V j 



r ■ 






' V ' 



» * • **■: i i 






/ u ■ 









f 



x v *.. # 

r >v. 



• ■ r i ' t t 






& 



j 



The correlation was plotted and compared with data a» tec nrcosurcs: 

90 psia (figure 6) and 2000 psia (figure ?). At 90 psia tie correlticn 
is outstanding and does an excellent job of predicting the data* At 
2000 psia the correlation is higher than the data, yet only by . 

•blrshak, et al: (ref; 13) 

(q/A) cr “ It §0,000 ( 1 + .0131 p )( 1 ♦ .0050b AT^Xl + .0365 V) 

Range of values: p = 25 - 85 psia, V = 5 - h5 ft/sec, D «* .21 - 
,lt6 in., whore AT S3 ^. is measured in deg F, and (q/A) cr in Bt)^/br-ft 2 . 
This correlation was compared with experimental data for ,2lt2 in. 
diameter tubes at 90 psia (figure 8) and .305 in. diameter tubes at 30 
psia (26). In both cases the correlation curve was conservative (lower 
than actual data). 

DuPont : (ref; 8 ) 

(q/A) cr - 595,000 (1 ♦ ,0ltV)(l + .0055 AT'satXl + .0055 ^ SCE ) 

Range of values; V ” 8.5 - i;0 ft/sec, = i+0 - 120 deg F, 

D * .3 - .88 in., p e 60 - 1000 psia, £T sc g * Subcoding Gradient * 
( ? exit “ "inlet) Per 2 foot length of tube. 

DuPont's correlation is similar to the one developed earlier by 
Eirsbak, et al (13)j a combination of linear factors of velocity and 
subcooling. In the DuPont tests a two foot length of tube was used. 
Subcooling Gradient (A r " 3C g), the temperature difference through the 
tube, was the last factor used to correlate experimental data, using 
the form of a previous correlation for annular geometry. 

Doroshchuck an d Iantsran:(ref ; 9) 

(qA) Cr - (q/A) 0 * A(p) [(.3l6/D)*33 - lj Qtu/hr-ft 2 

where: (q/A) 0 is determined from a plot of (q/A) cr vs. £h sat for 

8 mm (.316 in.) tubes for the pressure in question. 

Range of values: D « .117 in. - .316 in., p * 50 - 170 atm * 738 - 

25 00 psia. A(p) is a tabulated function of pressure. 
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p - atm psla 



A(p) - T>tu/hr-ft 



50 738 

80 111+0 



2.55 x 10 6 

: ,ia. x 10 * 



loo 1U70 



2.31 X 10- 
2.13 x 10 b 



1U0 % 2060 

170 2500 



1.99 x irr 



Based cn the 8 ran (.316 in.) data of Dcroshehuck and Lantsnan, the 
correlation predicts critical heat flux for the other diameters with a 
maximum deviation, of 30£, for fixed values of the other variables. * ’ith 
other experimental data the variation is slightly greater. 

KIT Small Piano ter: 



15 x 10 6 lbm/hr-ft 2 , i/D * 5 - 60, Ahg at « 0 - 150 Btu/lbm 
This equation was developed by the authors for the low pressure 
data from KIT and fits -ho majority of the data to i 2 $%» 



(q/A) cr - 1.738 x 10 6 (U.U8 ♦ .01 Ah ca t)(6/0 o )* 3i ’ (D/D 0 ) m *& 

where: G Q » 10 6 lbn/hr-ft 2 
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RESULTS Or i^P KDfctfTAL STUDY 



Scope ol Experimental Study 

As mentioned, there is a greet deal of experimental data available 
in the literature. Although each investigator has attempted tc be 
somewhat systematic in his own method of attack, there are very few 
cases where data M> from one investigation can be compared with that from 
another. Generally there are substantial differences in the variables, 
and lacking a good correlating equation, comparison, except in a very 
general way, is impossible. The experimental phase cf ttis study was 
conceived with this situation in rrind. Ornatski and Kichigin (li, 5 
and 6) have conducted tests with tubes cf two millimeter diameter, and 
four mass flow rates, over a range of pressures from 11*7 psia tc 281*0 
psia. Part of the present study was to complete the low pressure end of 
the 2 mm map. Also, using the same size tube and flow rate, other 
parameters were varied. The results of this investigation, when combined 
with the data from Ornatski and Kichigin, present a fairly complete 
picture of CHF for tubes of two millimeter diameter. 

Experiments were conducted to yield the following sets of data: 

a. CHF at various subcooling for three mass flow 
rates (G *> 3.69 x 10 6 , 7.37 x 10 6 , ll*.76 x 10 6 
lbm/hr-ft - identical to those used by Ornatski 
and Kichigin) with D = 2mm (.079 in.) for exit 
pressures of 30 psia and 90 psia. 

b. CEF at various tube diameters, with diameters 
varying from 0.0l*7 in. to 0.2l*2 in., with mass 
flow rate held constant at 3.69 x 10^ lbm/hr-ft^, 
and at constant exit subcooling of 35 Btu/lbm* 



2h 



i. With L/D held constant au a v ..lue 
of hO, 

ii. With L held constant at a value 
of 3.15 inches 

c. Gift for various L/D, for the same, masA flow 
rale and exit subcocling as in part h. above 
with D » 0.079 in. L/D from 3.5 to 60. 

d. CHF at various subceolings for the same 
mass flow rate, with L/D cf ho and 

D = 0.2h2 in. This data was obtained in 
order to present a comparison with the 
corresponding D * 0.079 data, and as 
substantiation of the GIF verses diareter 
data. 

e. CKF, for the same flow rate, exit sub-cooling 
and L/D, D * 0.09h in. with changing entrance 
conditions. In addition to the calming 
length form us- d in the majority of the 
study, four entrance forms were investigated. 
In all cases, the entrance was placed 
immediately up-stream of the heated length. 
Entrance forms investigated were: 

i. Sharp Edge Entrance, 
ii. hS-degr: e Chamfer, 

iii. Rounded Entrance, radius of 

entrance equal to 2 diameters . 

iv. Reentrant, length of reentrant 
section equal to 10 diameters. 
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f. ' 'axirum Heat FI jk. in ^ nvor r aMcn 
of the IT'S ter 1 capabilities , using 
0.212 in. T.n. fdunr, with l/ r ' of 
7»7i> and the **■•'■• J mun available flow 
rate arm 1 exit suhrocling, tc determine 
t;fb raximum obtainable C T ' at this 
diameter. 

Resets 

A total of 117 CHF data points were obtain©' • in .be experimental 
phase of this study. Of these, twelve were obtained usinr the Beryllium 
Copper tubing. A detailed description of the problems encountered 
with this tubing is contained in Appendix D. 

Figures 3 and it show the Of data obtained with I? «* 0.079 in., 
for the three :nnss flow rates. These figures show the same general 
dependence of CKF on mas 1 ' flow and subcooling c the majority of the 
other experimental data. This lev) pressure CRF data does show the 
-ininum CHI in the subcooled region more clearly than does the data at 

higher pressures. It can be seen that, n minimum CHI is encountered at 

approximately 30 Btu/lbm subcooling, and that the subcccling for 
minimum CHF increases for increasing mass flow rate. CHI 7 becomes 
nearly ind spend ant of subcccling at the lowest flow rate, at 30 psia. 

The dependence of CHF iron tube diameter ir shown in figure 5. 
Investigations were conducted with L/D held constant and with L held 
constant. The crossing of the two curves merely shows the diameter at 

which the selected values of L and L/D coincide. The data obtained 

with L held constant shows a greater CHF for the larger diameters, 
reaching a maximum increment of ap roximately 1$% at the 0.2h2 in. 



diameter tube. The value of CHF can be seen to become less dependent 
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upon diameter for the larger tubes. Limits of the apparatus precluded 
investigation of diameters larger than 0.2h2 in. at the flow rate 
selecte' ! fcr this study. 

The comparison between CKt data at D 3 0.079 in. ani that at 
D * 0.2U2 in. can bo seen in if pure This data provides an evaluation 
of the validity «i the c ur/e of CHT- verses diameter (f ip ore b ) which 
was determined lor a particular value ox exit subcooling. It can be 
soon that proportional difference in CHF between the two Diameters 
is nearly constant over the range cf exit subcoolings considered. 

Investigation cf the effect cl changing length-diameter ratio 
upon CHF (figure 7) shows that CUT increases for decreasing L/l). 

This data shows considerable scatter at the lew values cf L/f, 

The test section construction procedure, discussed in detail in 
Appendix C, was such that very small fillets of sclder remained at 
each end ot the heated length. In these small sections, this limited 
the accuracy of determining heated length. It can be seen that CFF 
approaches a minimum value as L/D is increased, and that r HF is 
essentially independent of L/D lor values of L/P greater than liO. 

This is in agreement with the results of Eergles (l, 11). It should 
be noted that the value of CHF for L/D of 10 (the generally reported 
minimum length in the Russian papers) is seme 1 $% greater than the 
minimum value. 

Results of the investigation of Entrance Efiect are not plotted. 
Holding all flew variables constant and charging only the form ol the 
entrance to the heated length yielded the following data; 
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Entrance Form 



Heat Flux 



Sharp Edge 
bS° Chamfer 



20- D Calming Length 



Reentrant 



Rounded 



3.65 k 10 6 Btu/ft 2 br 
3.86 x 10 6 Btu/ft 2 hr 
3.62 x 10 6 3tu/ft 2 br 
3.56 x 10 6 Btu/ ft 2 hr 
3.82 x 10 6 Dtu/ft 2 hr 



The value for the sharp edge is a mean of ,wc test values, the ethers are 
single tests. The 20-B Calming Length was the standard form used for 
the remainder of the experimental study, the value listed is the result 
of several tests at the sane conditions. Complete data for those tests 
arc lis'ed under runs no. 33 through no. 37 in the listing c.f data. 
Appendix C. 

In the sfudy of maximum heat flux, suhccc lings on the order of 
200 Btu/lbrr, and flow rates on the order of 5 x 10 :; lbr./ft 2 hr were 
obtained with a test section pressure cf 8p psin. These conditions lead 
to a CKT of almost 6 x 10^ Etu/ft 2 hr. Exact data cn these tests are 



contained in the listing cl data, Appendix C, runs through numbers 
29, 30 and 31. 
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IV CrvCLiJ yr T pis' M r fl TT^Io 

Discussion of Experimental Pe3ults 

The data obtained in the experimental phes oi -oh is stocks hews 
the same general dene ndance of C HF upon the ether variables as h«s 
been reported in similar studies. This data is perhaps a little more 
well ordered than that normally reported. This, of coarse, is 
heartening to the investigators, and is believed tc be the result of 
careful experimental procedures and constant cognizance oi the flew 
stability problem. It is felt that instrument and reaping accuracies 
would combine tc insure an experimental accuracy of five per-cent or 
better. 

The curves of CHF verses exit subcocling terminate at the maximum 
subcooling obtained. Thrs maximum is dependent upen the test section 
inlet temperature and pressure. The curves show the maximum subcooling 
available consistent with the capacity of the installed heat exchanger 
and the temperature of the city cooling water (minimum inlet temperature 
* h 7° F). The additional l£ or sc Btu/lbm of subcocling which could have 
been obtained by the installation of a complex test section pre-cooler 
did net seem tc justify the additional complication. 

The curves of CHF verses exit subcocling show a minimum value of 
CHF, signifying a change of flow regime. This minimum occurs in the 
subcooled region mainly because of the manner of determining conditions 
at the point of CHF. Fluid renditions were determined using standard 
thermodynamic var ables, which are valid for equalibrium states, but 
which dc net indicate the physical conditions with inherently non- 
equilibrium subcooled, forced flow boiling. What is needed is a better 
means of describing the physical state, one which would include some 
allowance for the non-equilibrium conditions as they actually exist. 
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The use of a vcid-fra 1 " oion variable heu x. wUg i o _ : : rue, 

(35) » but this is verv difficult uo G’/al x e>v .i -..pc if at v a ,-u.rh 
the sene variables as C’lr . The s? cf subetol? ng, as Jeter-: no- bv 
a difference in enthalpy or tcroe-sture, is ccnveri r . , r>u, hut 'uv„- re 
accompanied by a realization of its inherent inanility srr" : be 

the actual physical phenomena . 

At the very high values 0 . exit subrocling , bh c : e r , T 
seems to become less steep. That is tc say that the rate c change of 
CKF with subcooling decreases with ire-easing subccoling. This i very 
slight effect, not at all pronounced at low 37s ter pressures. The 
high pressure data cf horoshchuek and Lartsmn (10) and Ornriski and 
Kichigin (5, 6) shows this effect more clearly. A possible explanation 
for this lies in analysis of the effect cf subcrolinp uoon the ovors.ll 
boiling mechanism. It appears that once high values of subcoeling 
are reached', and a large recondensation potential is established, 
further ircreases in subcooling have less effect. 

The increase of CKF with decreasing tube diameter is evident in 
figure 5. It has been suggested that void -fraction increases with 
decreasing diameter. This is duo to the manner in which void-frscticn 
is defined, rather than peine physical effect. The increase in CPF 
is a direct result of the requirement for continuity cf mass flow. The 
size of a vapor bubble is dependent mainly upon pressure. In small tubes, 
the bubble occupies a proportionately larger part of the cross section, 
requiring that the liquid phase assume a higher velocity tc maintain 
continuity. It has been repeatedly shown, and can be easily explained, 
that CHP increases with increasing velocity? thus the increase in C5IT 
with decreasing diameter. 

The comparison of CKF with length-diameter ratio (figure 7) > 
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shows a leveliy c if a a value <. ? I,/" c <. ’-'royir r,e._ r i ' . ' 

a reasonable mi driun valu tl. t can * -mg 1 ne in exre ~irv m 1 Y 'L, 
where one is atte-pting to insure +h* 4 .ha efbct.: t; i.l * • L 

variables a separated. 

An interesting effect oi changing L/ r car. be s: on in 1. :,A f 

CHF verses diameter (figure 5 ) The curve for constant tube 1 oh 
re lects the effect of simultaneously changing -’iamet-.-r and 1 .n ,h- 
diameter ratio. At the larger diareter, the data far -t ns rent tube 
length reflects a substantial decrease in L/C, thus the greater CHF. 
Unfortunately, the Units of the apparatus were such that larger 
diameters could not be investigated at the mass flew rate selected for 
this study. It is altogether possible that the curve for constant 
length would show an increase in CKF for larger diameters; i.e. for even 
greater reductions in I./F. 

Comparison of CKF data for two tube diameters (figure 6), shows 
that diameter is one of the mere important variables. 7ithin the 
experimental accuracy of this study, one can conclude that the variation 
of CKF with diameter is nearly independent of subcooling, at least in the 
low pressure region. This requires investigation at higher system 
pressure, where the effect cf pressure beccres more pronounced, and 
where higher levels of subcoclinp are available. 

The data from the entrance effect study shows that entrance goonetry is 
unimportant in determining CKF, Here, as throughout this problem, there 
is a high degree of interaction between the variables. The investigation 
of changing entrance form in this study was done with a large length — 
diameter ratio to minimize the L/P effect. As a result, variations in 
flow due to the change in entrance form were likely damped out at the 
point of tube failure. Decreasing the heated length would have introduced 
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the compli'ctics. of the T/ T ' efJ'Tt. t in safe «,c r.n u;.e that c’jrnping 
eitr-.vce ’ C" bar no ef" -ct fcr i.i.h values of L/b. The scnevhet lower 
HT for ■■) j r- eat rant lurm re f reflect fit instability duo tc compressibility 
3 a the non-flow spare around the roan traat ser ticn, but 3s were likely due 
to expo ir.ental scatter. 

T) a rsximnn^h.eat flue data 'is cu interest onlv fcr future investications 
us5n t the particular apparatus employed in this study. This data has no 
value 3n a alysis of the overall HO problem because of the limitations 
imposed unen pressure and mass flew rate by the apparatus. Tests wore 
conducted to determine the maximum CHF that could be obtained with 'working 
size tube diameter (0»2u2 in.) an ; this apparatus. 
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H°< ,. n 'Tr’i ' * . ; . *; r»G j ' t t < 

The a c Hewing ''ornrcenie t n . ih various rar el; Hen*, a ■ pr 
srhenes are base-3 u or close examination err ccnparH • ith wq. 
data. The elect cf tho c< ntrollin poraretors the . . -t alluded 
described in detail in the section in Farmers T nil serein Cr'< ti< 
Feat Flux (Chapter III). 
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3oll: (rof; 21) 



The use ol this cor L ■ion if reswir wJ L nr ares Sj. <P'0 ;.■ ia. 

There is no restriction cn fh^ ,ujc diaoev.r, I'ut fr other reculvs, 

it has been shown CM has a strong i depend vie.? o r ’.iareterj 

therefore, this correlation should be ccnsi ?ered useful for laree 
diameters (0.25 "‘in, or lamer) only. 'Ithcuph the ac. uracy is to 
30 ^ in the case illustrated (figure 9), the cc^^lo tiers is relatively 
simple and c.-sy to apoly, and is generally n nservative in its prediction. 
Jens c- Lottes : (ref : 17 ) 

This correlation is of a power function fern, a»d wcult indicate 
zero CHF at zero subcccling unless a realistic lower lim* t «as 
established. Since the correlation das: gned to have nc negative 
deviation, it is one of the most conservative considered r ts main 
virtue is it3 ease of application. 

Ryabov and Berzina: (ref.’ 36) 

The 2 mrr correlation is adequate for the dais at pressures above 
lb7 psia. It is not consistent with experiment for the lower ore-; sure 
data (30 or 99 psia). There is 3 minimum subeocling, approximately 
18 Btu/lbm below which the correlation fails due to the notbe'-aticsl 
form. 

The correlation for D ^ Jtarcr is adequate for larger tubes, 
(approximately ^in. ), but for smaller tubes it suffers from lack of 
consideration of a diameter e. feet. Ryabov and Berzina considered the 
diameter effect in the derivation of this correlation by changing the 
leading coefficient, but they intimated that the diameter effect is 
negligible for larger diameters. This correlation could be improved 
by substituting a function of diameter for the leading constant. 
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Zenkevi ch; (re f : fh) 

IV ••luation cl viscosity A r hi. re- 1" ~r r „ d to*"- nature 

pro : uces a. curve sen what like tV? experiment d- t- ■> • is roc need 
if th" 1 ? ccr-cloui co is tc- be used. 

The absence of a tube dianeter factor limits ti.e uso.uVcs- of this 
correlation, Tfte lover limit, of 0.16 in. is toe small tc exclude all 
diameter effect. It is rocc-?- : r«rdoc that the use ct this correlation 
be limited to tubes of 0.25 in. in diameter or la" er, here cho 
dependence of CPI u on diameter is quite s r ai.l. 

Zenkevich r- Zubbot~ n : (refg lb) 

This cor 'elation, like so "any ethers, ner loots the effect of tube 
diameter. The lower limit cf 0.16 in. is too small to neglect this 
effect entirely. For the larger diameter (.316 in.) data; this 
correlation is excellent. The lower limitation on aubccoling is sound, 
due tc the math-enstieal fem of the equation. The upper limit has 
no phrr'ca 1 validity and seems t . be evidence cf caution. In spite 
cf the diameter ccnsiderao'cn, the mathematical expression is 
convenient, and tie results for larger tubes arc quite rood. 

Gam, bill: (ref: 30 ) 

The curve for wall tempera ure is suspect, and available nata of 
wall superheat at burnout ( ce s net bear out the graph presented. There 
is too much latitude in choice of constants and exponents to bo used 
in the correlation. Presumably, if the data is reasonably veil ordered, 
the fit of the cor"- laticn ran be improved, by adjusting the values of 
K, K* , ra, and n. There are no easilv applied guidelines for the 
adjustment cf these cons tarts included its the presentation. This 
correlations does include a diameter ter", but it is only in the forced 
convection part, and then tc only the (n - l) power, where r» is me of 
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Gambill’s "variable" cons ants, generally on tho order of ^.3. This is 
one ol the few correlations which includes snail tube diameters, yet 
clearly the correlated effect of diameter is loss than that observed in 
the data of I'ergles (l), Poroshchuck and Lantsman (9), and figure In 
The correlation is no better than that of "irshak (13) at la* pressures 
and is considerably clumsier to handle. At high nressuros (figure 9) 
this correlation is no better or worse than many ethers. 

Be math : (ref: 16) 

For a given pressure, tube diameter and flow rate, this correlation 
is linear with subcooling. This is a faijfcly good first approximation 
and typical of many correlations. This correlation is rather dif ficult 
to apply, requiring auxiliary calculations tc determine heat transfer 
coefficients and wall temperatures. Of the seven high pressure 
correlations plotted on figure 9, this is one cf the best, predicting 
CKF within for the selected conditions. 

Bettis > ( ref : 17) 

The use of this correlation is restricted to pressures near 
2000 psia. There is no restriction on a minimum diameter; therefore, 
the correlation will yield low predicted values of CKF for diameters 
less than .2!? in. 

Griffith: (ref: 13, 17) 

This correlation was ogiginally determined for several fluids 
including water and it is applicable throughout the subcooled boiling 
region and into net steam generation. Generally, over a wide range 
of variables, this ccrrelatj.cn is the most consistently accurate one 
available. Its disadvantage lies in its extremely cumbersome and 
unwieldy form. Computation can be simplified somewhat by using the 
curves presented in WAPD-188 (17) for several cf the recurring variable groups. 
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? arshak, et al: (ref; 13) 

The for:- cf the ccrrel at ? is Linear X X rccli ®r u 
exhibits the approximate tendcnr” • f t : '- er<r>f ri itrtal d*L» . ‘it* in 
the range v.f diameter.-, icr vl i < h «ho cornel.- thn : s , 'he 

depet. denre cl OUT cn -'is ns ter is neglects' 1 . The Icllo-inp g -serai 
cements are appropriate : 

1) The correlation Ives a cc nservat ive In 'L-ation 
of nit, but errs by 25% - sO'" 

2) There is no consideration of the variaticr of 
CRF with diameter 

3^ The correlation is safe, ouick and e^sy to amir 
^jPont: (reft 8) 

This correlation was designed fr~ correlating r, u T ‘o' data for 
tubes of a specific length. How to evaluate Tub^colin-- f-rrPient for 
ether tube lengths is unresolved, ^ubccoling gradient is not an 
independent variable in the equation, being tied to Heat I lux by the 
First Law of Thermodynamics, lia meter, ~nd flow rate. Ha diameter 
effect is noted between 0.30 - 0.88 in. This is a reasonable 
approximation. Variation in pressure between 60 -1000 psi3 is 
considerable and should be taken into account. Since the equation 
was derived from a snail sample cf data and for tubes of only one 
length, its use is net rec emended. 

Doroshchuck & Lantsnan: (ref; 9) 

This is one of the few studies which has been devoted exclusively 
to the diameter efiect problem. The general trend is correct, but 
comparing this correlating equation with experimental data - even that 
frer the paper in which this cor elst'cn is presented, shows an error of 
up to 30$. Part of the difficulty is due tc the interaction of the variables. 
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Covpari ng t’ is rrhr"3 •. itb c her e < re ' d-tn, nt«ably that 
of c rratskl rnc 1 "<ich'"in (h, f ~ , ) -rd ’■ r< los (l, t} ,) "tc shored 

the ccr-’f let 4 cn tc exceed a lOfi e w or. ; ; 's zo.iu vhr a cf 

A |j(D/D 0 )^*" - lj fit the Orootski 2 nd Kirhirin data (i;, f and 6) 
quite well, but this failed in ccrmsriscn with dot". frt~ otter s . dies. 

It is net r^ccmmonded that this corbel? on be ir-~d tc ’^d'et 
Ph'F j hr never , it will serve tc arnly a di are consideration t c ether 
correlations, within rather bread Units, 

'IT 1ml I r lameter 

Th's equation was develooe ! not sc irivh a ore die icq of r Ht , bat 
tc provide a fairly r i noise mathematical expression for the leu pressure 
t 7 !? J ata. ? Tc pronouns tern is i- laded as the variation of HF with 
pressure is neither well ov<3eve> A nor of la-rr-e magnitude. The change 
in TFT -.r th pressure ever this range is "rich less than the o~r or of 
this equation due to the assumption cf r?J as a linear function of 
subccrlirr. nnr is, in fart, not linear with sub^cclinp, this form 
has been assumed merely for convenience. Tt is recommended that the 
use of *his equation be limited tr deter -ination of v< neral levels, 
and that the actual data be used in any attempts at accurate prediction of 
CHF. 
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S* r LL \ v Tf'R r?IT' r " T - H 'ft. ' ’ Lj ” j_ r T -. r " h£.r 

PKCC njRC HP ugv 

Figures 3, it, 10, 11 and 12 provide a rraohical present a tier of 
CIIF for pressures from. 30 psia to 2oitO psia ever ? range of mass flow 
rates and exit s ^coolings for tubes of two mn diameter. These curves 
in conjunction with the data of CHf dependence upon tube diameter 
provide a means of graphically determining CPF for tubes uo tc one- 
fourth inch diameter, flow rates from 2 x 10^ lb m /ft" hr. to 22 lb r(1 /ft^ hr, 
at any pressure and exit subcorling. A procedure of linear and graphical 
interpolation is recommended for predicting r HF for any reasonable 
small diameter application. In cdses where linear interpolation is 
recommended , auxiliary curves, have shown this technique to be 
satisfactory. The accuracy of this graphical prediction technique 
will exceed that of the few correlating equations applicable to small 
tubes, and in the majority cf the cases tested, was within Vo% cf the 
experimental data considered. 

The following procedure is recommended for predicting CHF : 

1. Select desired system pressure, mass flow 
rate, tube diameter and length. 

2. Locate the plots for pressures bracketing 
the desired pressure. 

3. At the selected exit subcocling, read 
CHF for the three or four mass flow rates 
presented, for each pressure. 

it. Prepare interpolation curves with coordinates 
of CHF and Mass Flow Rate, and determine CHF 
corresponding to the desired flow rate at each pressure. 
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5. Using the tuc values thu- ortainco , linear 
interpolation between th' uroscures un-’er 
rc ns id oration will f i ir e a value cf C'F for 
the c.esired pressure, flow rate and 
subcooling, for a diaT.etcr of two ran. 

6. Refer to figure 5, and using the curve for 
constant L/P. determine the value of CPF 
for the ’iameter in question. The ratio 
of this last value, tc that at two ram 

(a value of CHF of b. 1? x 10 6 Btu/ft 2 hr, 
at 2 mm. diameter will provide satisfac ,crv 
results) multiplied by the previously 
determined value of CPF, will yield a 
predicted value of CKF at the desired 
conditions cf pressure, flow rate, subcocling 
and tube diameter. 

7. In cases where the desired L/D is les- than 
20, the same sort of ratio correction should 
be performed, using figure 7. For values of 
I/D of 20 or greater, the error in CHF vail 
be too small tc warrant consideration. 

For a rough, quick estimate cf CHF, visual interpolation for mass 
flow rate on the plot for pressure nearest the desired pressure can be 
easily accomplished. The diameter correction should be included in this 
abbreviated procedure. This same sort of approach is useful in 
estimating the effect of gross change in the important variables. 
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Figure (l) SCHEMATIC LAYOUT OF EXPERIMENTAL FACILITY 
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Figure (3) - Dependence of Critical Heat Flux on Subcocling 
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Figure (5>) - Dependence of Critical Heat Flux on Tube Diameter 
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Ficuro (6) - Popenderre of Critical Foot Flux on ^ubecclinc 
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icuro (7) - ^crtfcnra of Critical Teat Flux on Nratod Length 
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Flfjuro (0) - Comparison of Correlations an^ I5cocr:lnontnl 
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Figure (10 ) - ^opendenco of Critical Hont Flic: on ^ubnooling 
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Figure (12) - fepcndonce cf Critical Font I'luftc on ^ubcocling 
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Exp; rinonta on was rcnductr 5 on ex-.r- r ~ : r "at. is w' ch had 
been constructed i cc a previous ssuiv (7.^ •’ 'Jr'*- rat- been -is d 

for several other belli* r and rJT investi** ‘ icn . "'nly miner 
modification was ne essary for this investigation. The apparatus 
is shewn schematically in figure 1. The fluid loop consists of a 
closed rain circulation loop, with a by-pass parallel test section 
line. All components are ol corrosion resistant materials. The 
system contains a main circulation pump, acumulatcr, the test section 
line with flowmeter, preheater, test section with instrumentation, 
main loop by-pass line and heat exchanger with city water coding. 
Auxiliary to this are a fill pump, supply tank, degassing tank am’ a 
continuously operating de*in©ral*i*er. Fewer ir supplied to the test 
section by a pair of a.c, motor - d.c. generator sets. 

The main supply pump is a two- stage turbine type pump, driven by 
a 3 hp induction motor. The purp i* capable of delivering 3.6 gallons 
per minute at 260 psi. A bladder type accumulator located, at ihe 
exit side of the purp serves uo 'amp out pressure fluctuations. Plug 
valves at each end of the by-pass lino provide control for the test 
section pressure. Test section flow rate is determined by a Fisher 
Porter flowmeter and controlled with a Hoke needle valve. This needle 
valve also provides the 13' 0 psi up-stream pressure drop necessary to 
insure flow stability. 

The mctcr-gcne rrt or sets consist of Uij.0 v. , 3 phare synchronous 
motors and 36 kilowatt d.c. generators, connected in series. Each 
generator is nominally rated at 12 volt and 3000 amps. A total power 
of approximately 60 kilowatt is available at the test, section. 
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Instrumentation is available To cnitcring pressure levels, 
tempera cures, test section flow rate and test section current and -dtage. 
Test section exit pressures are read cn bourdon type pressure gauges 
with specified accuracy of 0.£ psi or better. The gauges are calibrated 
periodically with a dead weight tester. Test section flew rate is 
determined with the Flowrater ’etersj the range of flew rate is varied 
by changing tubes an. floats. Meters are calibrated in place and 
calibration curves are available. Test section power is obtained from 
the test section voltage and current. Voltage across the test section 
is read directly with a ,-Jestcn variable scale, calibrated d.c. voltmeter. 
Current is read across a National Bureau of Standards shunt, with a 
specified calibration of 60.17 amps/r.iv. 

All temperatures are determined by thermocouples manufactured 
from 30-gauge, duplex, Ccpper-Ccnstantan thermocouple wire. Test 
section inlet and exit temperatures are measured on thermocouples 
inserted into the flew with Conax fittings. Output voltages of the 
thermocouples and current shunt are displayed upon a Minneapolis- 
Honaeywoll Brown, 0-26 mv., single channel recorder. 

A much more detailed description of the apparatus is contained in 
ref. 1 $ however, certain modifications have been made. The more important 
of these include: 

a. The braided copper flexible power connection 
and the rigid up-stream jumper from the 
test panel to test section bus connection 
were replaced during this study by Kackwcrth 
Reese flexible Dry Jumpers. This requires 
that the test section be supported. Support 
is pervided by a laboratory support and clamp, 
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•with elect rical irs\u->.„ i. 

b. The braided copp~r peer supply cables have 
been •'•eplacrr' .-ith v arkr'ortl 'cose , water 
cooled cabins. Theee cable o an. cooled by 

city water tapped off c •' the ho? o;cha: per 

»> 

coolant supply line, discharged to a drain. 

c. The preheater has been ~odified by the 
addition of three $ kilowatt imrersion tenters, 
with '••operate controls. These heaters are in 

serides with the original 6 kilowatt variable 
control heater and provide any desired 
preheating from 0 to PI kilowatt, with 
steady state operation. 

d. The stcan heater in the supply tank has 
been removed as it wan found that 
sufficient heating and degassing could be 
obtained from the degas tank heaters alone. 

e. The thermocouple selection switch now 
consists of two Leeds end North rap, 12 
position selector switches, in s ries. 
Thermocouple wires ate connected through 
screw tvpe terminal boards, thus reducing 
set-up tine. This modification was trade 
in order tc incre- se the number of 
thermocouples whir’.- could be patched 

to the recorder. Bcsictarce tests showed 
no appreciable change in resistance as a 
result of the added switch. In this study, 
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only th ''in? ■•■ svi^ ~ » s u. zoo., 

i. When using -.ho vor t Short ter: sections, 
nr there cl tvpe r the test 

section recistmeo ra:- *;v lc a:- tv 
’(Bke the rc tor- goners ter 
highly erratic. To alleviate this 
prcbler, a series resistance vies 
fabricated. This resist "-nre consisted 
of a 10 inch length of 0.2ii2 inrh type 
30’; stainless steel tubi ng, with bus 
ccnoeotionr, place 4 :lectr J c->Hv in 
series with the test oc-c+ic'- and oocled 
b/ ofty water. This added resistance 
caused an additional voltage d**o 0 ? 
ap r c coin te 1 y *> volts acres? the poor 
supply, which im roved the r,-g. set 
ccrtocl considerably. Power to the teat 
section was deter^i^ed by readier voltage 
drop across the test section cr'ly, and by 
reading current through the series 
conhi nati cn. 

g. In the tests whe^e verr high power 

applications were ■’-equired tc i*carh PHF, 
the tost section exit bus connection was 
found to heat up excessively, causing 
sacking and in '•nr ca~° a nr-all fire, not 
tc nenticn erect i c'-^nvenience when chanyin 
test sections. Tn v'ew of this, a bus 
connection coolant scherzo was devised. A 
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• -inch copper n: _e r two el’- us van 

silver brazed tc the brass bus connection. 

Cit ,r water, fror the ’"■oner lend cooling 

system, was diverted throuf h ohls rdnplu , 

drawing herv from the bus connection. Heat 
£ 

balance checks when operating, comparing 
power input to fluid temperature rise through 
the test section showe 5 that this cooling 
device had no e: feet on the test data. No 
further overheating problems were encountered 
once this cooler was installed, and it was 
used whenever there was a require ent fer 
high power input. 

Test Section 

Because each test destroyed the test section, construction was kept 
as simple as possible. A test section was constructed bv facing a j/h inch 
brass pipe cap and drilling a hole in the center of it. A section of 
the test tubing -was silver soldered into this hole. At the ether end 
of the desired heated length, a 5>/l6 inch brass bushing, about 5/8 inch 
long, drilled to the O.D of the tubing was silver soldered to the 
tubing. Heated length was measured from fillet bo fillet. A calming 
length cf approximately 20 diameters extended beyond the up-stream 
bushing. 

Sections were installed in the apparatus by threading the pipe cap, 
the down-stream end, to the down-stream plenum chamber which contained the 
exit pressure tap and exit thermocouple. The up-stream end was installed 
by inserting the calming length into a Conax fitting sealed with a 
neoprene washer. Brass bus connections were installed on the bushing and 
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pipe ca; rcm v 'l tir Vic ( s c + 

In the st ad" c .uxi »„ " • .< -• *c . --hir.- u.> av< 

diameter cl tlx. C,2i.Z in. :. v i . o' , r c -••r t. cM/u- 

adequate impcndai.ee ncc-c'd a 1 .Mr in ’ .c. i; noted in the 

15 stinf or 7>ata. (.‘'nponiix ", ru/o ?, 3C : o 31). 

% 

Be rvll iu~ - C cope r Tuhiiy , 



At the be? irdr. t *J the env "inentc lien, 
was erde^e with i.usidc- diameter ef P.O' 7 9 i 
that this tuhiri; -■ ai ho or - ic» t* • h t 
steel which hat hoc as i in r-: - v.c ..t ' -'ic 
vers nsado up nnt teste rer« run. •rt’xar.t 

was encountered usin'* the Ec-fu u v v . °r. ’ 

problems are discussed below . 



i i :rv cj 2? r ■ lop, -le- r u 
nrrcc ( c m'. It was felt 
’ f Ivr. tha rpo 3°L stainles 
. isrt sections 

el..', c; c:"Hr’*aKLe diffi ultv 
■ x cf the : r - ’ important 



About the only iabricavl.cn advents 3 thn 1 the r e- r u tuning had wa 
that it could be cut _nto .desire 4 lenr'ir. -'l’;. unfortunately, the 
alley was not s fC j xo no t*-.. p rat-ares an excess of 1100 degrees . 

T f test sections were constructs' 4 by silver soldering, th tu.irc near 
the soldered cc.. notions underwent a phase change a:. ' this was a- conpanied 
by a change in electrical resistance of alr.cst 1*0, . Since this cliange 
in resistance oc erred only in the region where heat was applied for 
soldering, the lcoal c-Mvre in resist" nr o caused an uneven arcia? 
distribution cf halt five. t 'at-' x cbt-nlrod b~ ws*»!«»uri n~ tctcl to-t 
section voltage drop and current did net indicate the actual HI . 

Test sections were constructed ucinr soft solder, but the hirh 
exit temp ratures er.c curt ore' 4 c’-iri^g o oration on ised the soft solder 
to puddle gust before reaching Cff . «Sc lealcs developed, bat this was 
felt to be less than sat is far term. higher temperature (bO-lo) soft 
solder was tried, with no better results. 
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The lower r. sistance Be-Cu tfoi-%. reused m.g. set rcntrcl problems. 
The series resistance previously mentioned alleviated this nroblnn, but 
tended to complicate the procedure. 

HE- data was obtained using t 1 is he-Cu tubing. 'ubsequent comparison 
with similar data obtained with type 30b stainless steel tubing shewed 
that the Be-Cu data was inconsistent, for either the soft soldered or 
silver soldered sections. 

A corrosion check was made with the Be-Cu tubing. After one hour 
at a heat flux of about one-third of CHI , the inside of the tube showed 
distinct blackening, using derassed and demineralized water. 

In view cf the various difficulties, it was concluded that the 
Be-Cu tubing was unsatisfactory for experiments of this type, and thick 
wall, type 30b stainless steel tubing, with inside diameter cf 0.079 
inches (2 mm) was obtained. 

The stainless steel tubing could be silver soldered satisfactorily 
and the data obtained from it showed considerably more consistency. 

The data presented as the result cf this investigation was obtained using 
type 30b stainless stoel, except where noted to be Be-Cu (runs f 2 through 
013) t and except for the D = 0.180 inch tubing which was type ’’A" Nickel. 
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_ >_ i : .C “f Of 

lea/cs , 1* .-r ^ ^ig j. *_ _ * -i.3 l i j. ■->. .' »/or . 

?et-ccck£. . . ir jlo«. -i .n. lei .-it- .« w u • ; into throughout 

the a opera t us, an*. .-ir v: r hlec fr.r il & v ste . , both under the ../at 5 c 
head oh the • toga. 3 c nk, • <».' vico tie p-irn ir- ul; t.inf twcr through 
the loop, ^isolved air otlu .2 t ;r* .as ah*n rvffwta . rc" &h« test 
watii* electrically ho ti.r.r „•. ..<fcoinf oa>ik. As -sent Iciieu ir, 
ref. la.’d 2, about ov.c hal. ho. «. . bo c ;& «ai> •ecoeaar.y tv insure 

s at is f artery air rcrcval. 

C-xu'ti the do- acsinp ■ race dure was fo~pjjatert, an. he -water i tnc 
1. op ccolo- hf the boat exchanger, oh- dteirsd teat section flow r?i.e 
and .nr' -iru'-fs .-ore o. tsblioheu. 

Since it ».sa cesire . to pres nt daw. cr th ad a particular exit 
subccoling, it wae necessary to d'-tcriae and establish a desirec. test 
section exit „erpur* ourc. ills r-er rcrpiirated ay the iga that e. it 
torpor.- tore is tependent u...cn the inlet Lance rssure, which eculd oe 
controlled by the- preheater, and upon the nest a. lux in trie test section. 

In order to obtain the desired exit temperature and pressure at the 
point of r n , it war necessary tc i nc rea c the ower to the test section 
very slowly, and tc continually -.d ju. o the tecs section pro sure *nd 
the preheater controllnp the inlet temperature. In scute ewe, the 
21 kilowatts of preheat available was net sufiicient, in whi' h cas^ the 
flow of coolin' waier to tie hca« excrau; ?r was reduced. This procedure 
was diff ir ult in that cue exit tenner." -ure showed a long (about 2d minutes) 
response tire to adjustments in th- flew o- ity water ccoiinp the l low 
threw h tho he a./ exchanger. 
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Often the rocess of obtain- n the lesxrod exit subccolirtf .is c 
very delicate one. This was particular'/ true when ‘he rurvtf* of 'Hr 
versus subcooling was approximately parallel to the eurve cf. heat flux 
versus subcooling as determined ’ey the first Law of Therrodynami.es. 

In the large majority of the tes~ runs, it was possible tv. obtai. 
the desired valus*. of exit subeooling within 5$ at the point of CHF. 

In these few cases where exit subcooling varied mere than $% from the 
desired value, an auxiliary plot ci CHF versus subcocling was made, to 
determine the appropriate value, or to determine that the variation of 
CHF with subeooling wrs not si, nificant in the range cf subcooling in 
question. 

For the purpose of this study, CHF was defined as the heat flux at 
tube failure. While the heat flux thus determined may exceed that at 
the departure from nucleate boiling (DNB), heat flux sc defined should 
coincide closely with the heat flux at the point of transition from 
nucleate to film boiling. No burnout protection device was used in 
this study, as it was felt that less time and materials were involved 
in the fabrication of a new test section for each run than in the 
construction cf such a device. Furthermore, the^e rm be no doubt 
that a genuine maximum heat flux has been established when the tube 
ruptures . 

All experiments were conducted with the test section in the 
horizontal position. System water was degassed only periodically. 

Care was taken to insure that no air was trapped in the system when 
replacing sections and extensive air bleeding was done before commencing 
a new run. The inside surface of the test section was given a thorough 
cleaning with acetone prior to installation. No additives were 
introduced into the water. 
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In almost all cases, .-he tu e ruptured witt * n a ree c’i me. tors 
of the exit connection. Thus conditions in the exit pie ur chanb )r ■'•cro 
considered to be representative o_ those at ohe ocint c< failure. A 
recent study by Lopina (3) has shovr that the pre sure within tt" plenum 
is virtually identical with that at the down-stream end of the h nted 
section throughout the subccoled region. The test section and plenum 
were wrapped with fiberglass insulation to greater than the critical radius 
to minimise heat loss. 

Because the recorder used presents only one channel, it was 
necessary to estimate certain quantities at the point of CHF; however, 
near failure, very small increments of newer we^e an oiled , and all 
quantities were checked frequently. Bata not obtained right at failure 
could be readily inferred from that recorded just prior to reaching 
CHF. All changes in conditions were made very gradually during a run, 
and once the procedure was well established, it was possible to obtain 
a set cf data, e.g. one data point, in approximately forty-five minutes. 
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APPENDIX ? - SAMPLE DATA SHEET 



D 0. 0*1 In. P $ 0 _ _ psia Run # 45 

L -5. lb> in. L/D 40 Date z/z s 

G 5.69 X \ o 6 lbm /f t 2 hr = w = 125.S lora / hr 

= *74 $ on Flov/rater 



Flow % 


Tin ®v , 


Tout rav , 


P psig 


, E v 


I mv ,Variac ^ 


74 


Cc 71 '«*P> 

0.42 


£ 14-9 *P 0 

2.70 


ns 


695 ” 


3.52 


o 


s' 




5 . 04 


*•» 


7. SO 


9.06 




\S 




3.42 




3 .oo 


9.^6 






is' 


3.72 




3.40 


, 992 


- 






4.20 


76 


1 9- os 


10.60 






t ^ 
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Keat Balance 

wAT = 12.5. S *92, - 12500 

E«i( 6 o.i 7 )( 3 .Ai ) = < 2 . 95 * 3 . 62 * 60 . 17 * 5.41 - 12 ISO 



CKF = B • I( 60 .1 7 ) ( 3 ♦ A1 ) 
ttD-L/t 44 



CKF = 1O»^2.CS»(60.17)(3.41 ) 

TTvO-Ol^ZiS / 1 44 



h sat 290.3 
h bulk 223. Cg» 

hsat 4.1.7 Btu/ rD!n 

= Btu/ j. t 2 



Obtain Temperature from Thermocouple Conversion Chart 
h sat = h at p h bulk = h at T ou ^ 
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APPENDIX G 



DEVEI/)PHENT OF THE HIT SNAIL DIAMETER CORKEIATION 

This equation was developed by the authors for the low pressure 
MIT data and fits the majority of it to - 25%, 

(q/A) cr « 1.738 x 10 6 (U.U8 + .01 Ah^) (0/G o )*35 (D/D 0 )“*#^j-.12 

The correlation was developed in the following manner: The 30 

and 90 psia plots (figures 3 & 10 for 2 mm data (.079 in.) were linear- 
ized to determine the slope and the intercept of the part of the equation 
dependent on subcoolong. The flow rate dependence was determined by 
plotting on log-log paper, cross-plots of the 2 ran data at constant 
subcoolings . Also plotted were cross-plots of Bergles'(ll) and 
Vfessel's (26) data for constant subcoolings. A mean value of the slopes 
was used to determine the exponent of the flow rate term. The same 
technique was used for diameter and L/D dependence, using the data of 
figures 5 and 7* and the data of Bergles (11). The leading constant 
was determined by the zero subcocling intercepts and the values of 
the constants G 0 and No pressure term is included as the variation 
of CHF with pressure over this range is neither well ordered nor of 
large magnitude. 

The correlation was plotted over the KIT data and small adjustments 
were made in the constants and exponents to get a better all around fit 
and balance out the positive and negative deviations. 

A f urther refinement of the correlation would be to make the 
heat flux a weak 2nd or 3rd order dependence on subcooling, eliminating 
the linear approximation. A minimum could be introduced at about 
A^sat “ 20 to 30 Btu/lbm. 
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APPENDIX H - NOMENCLATURE 



a - Acceleration in g’s 

A - Heat Transfer Area, ft"" 

c - Specific teat, Btu/lbm-deg F 

D - Tube Dfajaeter, inches 

g - Acceleration of gravity, ft/ sec" 

gQ - 32.2 lb^-ft/lbm-sec 2 

h - I-feat Transfer Coefficient, Btu/hr-ft^-deg F 

h - Enthalpy, Btu/lbm 

J - Mechanical Equivalent of teat, 778 ft-lbj/Btu 
k - Themal Conductivity, Btu/hr-ft-deg F 

L - Tube Length, inches 

p - Pressure, psia 

Pr - Prandtl Nuidber 

q/A - teat Flux, Btu/hr-ft^ 

t - Temperature , deg F 

T - Tenperature, degrees Rankine 

AT 3at " ^ T sat - T b)j deg F, Subcooling 

- Specific Volume, ft-'/lbm 

- Velocity, ft/sec 

- Viscosity, ft'-/hr 

- Viscosity, Ibra/hr-ft 

- Density, Ibm/ft^ 

- Surface Tendon, lb^/ft 

- (h sa ^ - h^), Btu/lbm, Subcooling 



v 

V 

r 

<r 

Ah 
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Subscripts 



bo - At the point of burnout cr critical heat flux 

cr - At the point of burnout or critical heat flux 

b - Bulk or mixed moan conditions 

crit - Critical 

fg - Saturated Condition - Vapor Condition 

g - Vapor Properties 

1 - liquid Properties 

pb - Pool Boiling 

o - Saturated Conditions 

sat - Saturated Conditions 

v - Vapor Properties 

w - Hail 
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